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ABSTRACT. Selected results on estimating cosmological parameters from 
simulated weak lensing data with noise are presented. Numerical simulations 
of ray tracing through N-body simulations have been used to generate shear 
and convergence maps due to lensing by large-scale structure. Noise due to the 
intrinsic ellipticities of a finite number of galaxies is added. In this contribution 
we present our main results on estimation of the power spectrum and density 
parameter fl from weak lensing data on several degree sized fields. We also 
show that there are striking morphological differences in the weak lensing 
maps of clusters of galaxies formed in models with different values of fl. 



1 Introduction 

Weak lensing by large-scale struc- 
ture (LSS) shears the images of dis- 
tant galaxies. The first calculations 
of weak lensing by LSS (Blandford 
et al. 1991; Miralda-Escude 1991; 
Kaiser 1992), based on the pioneer- 
ing work of Gunn (1967), showed 
that lensing would induce coherent 
ellipticities of order 1% over regions 
of order one degree on the sky. Re- 
cently several authors have extended 
this work to probe semi-analytically 
the possibility of measuring the mass 
power spectrum and cosmological 
parameters from the second and 
third moments of the induced ellip- 
ticity or convergence (Bernardeau et 
al. 1997; Kaiser 1998; Stebbins 1996; 
Jain and Seljak 1997; Schneider et al. 
1997). 

Analytical calculations suggest 
that nonlinear evolution of the den- 
sity perturbations that provide the 
lensing effect can significantly al- 
ter the predicted signal. It is ex- 
pected to enhance the power spec- 
trum on small scales and makes 
the probability distribution function 
(pdf) of the ellipticity and conver- 
gence non- Gaussian. We have car- 
ried out numerical simulations of 
ray tracing through N-body simu- 
lation data to compute the fully 
nonlinear moments and pdf. Details 
of the method and results are pre- 
sented in a separate (Jain, Seljak & 
White 1999); here we summarize the 
method and present some highlights 
of the results. Other recent numeri- 



cal work includes Wambsganss, Gen 
& Ostriker (1998), Premadi, Martel 
&; Matzner (1998), van Waerbeke, 
Bernardeau & Mellier (1998), Bartel- 
mann et al (1998) and Gouchman, 
Barber & Thomas (1998) 

The dark matter distribution ob- 
tained from N-body simulations of 
different models of structure forma- 
tion is projected on to 2-dimensional 
planes lying between the observer 
and source galaxies. Typically we use 
galaxies at z ~ 1 with ~ 20 — 30 
planes. We propagate ~ 10^ light 
rays through these planes by com- 
puting the deflections due to the 
matter at every plane. Fast Fourier 
Transforms are used to compute gra- 
dients of the potential that provide 
the shear tensor at each plane. The 
outcome of the simulation is a map of 
the shear and convergence on square 
patches of side length 1 — 5°. Sev- 
eral realizations for each model are 
needed to compute reliable statistics 
on scales ranging from V to 1°. 

In the weak lensing regime, the 
magnification and induced ellipticity 
are given by linear combinations of 
the Jacobian matrix of the mapping 
from the source to the image plane. 
The Jacobian matrix is defined by 



where 60 i is the i—th component of 
the perturbation due to lensing of 
the angular position on the source 
plane, and 6j is the j—th compo- 
nent of the position on the image 
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Figure 1. The dimensionless power spectrum of hi. For the cosmological model 
indicated by flrn ctnd T in the panel, the power spectrum from ray tracing 
shown by the solid curves is compared with the linear (long- dashed) and non- 
linear analytical (short- dashed) predictions. The angular wavenumber I is 
given in inverse radians - the smallest I plotted corresponds to modes with 
wavelength of order L ~ 3°, where L is the side-length of the field. 



plane. The convergence is defined as 
K = -(^11 + <l>22)/2, while the two 
components of the shear are 71 = 
-(<^ii -<^22)/2 and 72 = -^12- The 
convergence k, can be reconstructed 
from the measured shear 71, 72, up 
to a constant which depends on the 
mean density in the survey area. If 
the survey is sufficiently large and 
there is little power on scales larger 
than the survey, this error can be ne- 
glected. 



2 Results of Simulations 

The N-body simulations used for 
the ray tracing are taken from four 
different cosmological simulations, 
the parameters of which are summa- 
rized in Table |^. The N-body simula- 
tions use an adaptive P^M method 
with 256^ particles, and were carried 
out using codes kindly made avail- 
able by the Virgo consortium (e.g. 
Jenkins et al. 1997). Coupled with 
ray tracing on a 2048^ grid these sim- 
ulations provide us with a small scale 
resolution down to ^ 0.5^, well into 
the nonlinear regime for weak lens- 
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Table 1. Summary of the parameters used for the N-body simulations, h is 
the Hubble constant in units of 100 km s"-*^ Mpc"-*^, F is the shape parameter 
of the power spectrum, and the other parameters have their conventional 
meaning. 



Model 




Qa 


h 




r 


SCDM 


1.0 


0.0 


0.5 


0.60 


0.50 


rCDM 


1.0 


0.0 


0.5 


0.60 


0.21 


ACDM 


0.3 


0.7 


0.7 


0.90 


0.21 


OCDM 


0.3 


0.0 


0.7 


0.85 


0.21 



ing. 

The power spectrum of hi mea- 
sured from the simulations is com- 
pared with the analytical predictions 
Jain &; Seljak (1997) in figure 
We will focus in this contribu- 
tion on the extraction of parame- 
ters from simulated noisy data. Fig- 
ure shows the reconstructed con- 
vergence field from noisy ellipticity 
data on a single 3° field. The galaxies 
were randomly distributed and as- 
signed intrinsic ellipticites with each 
component drawn randomly from a 
Gaussian with rms=0.4. The recon- 
structed hi is compared with the field 
without noise, as well as with a pure 
noise field, to show the significance 
of reconstructed features. The field 
is smoothed on the scale at which 
the noise and signal should be com- 
parable. The figure shows that with 
observational parameters feasible on 
large telescopes, a field a few degrees 
on a side can allow one to reconstruct 
the large-scale features dominated by 
groups and clusters of galaxies. Sta- 
tistically one can do much better. 

Figure ^ shows the power spec- 
trum of hi measured from simulated 
noisy data. On scales larger than 
10' the signal dominates the noise. 
The scales on which the power can 
be measured are dominated by the 
density power spectrum on scales of 
about 1-10 h-^Mpc at z ~ 0.3. 
Thus weak lensing surveys a few de- 
grees on a side will be sensitive to 
the dark matter power spectrum on 
these scales. 

Figure ^ shows one attempt at es- 
timating the density parameter Qm 
from simulated data. The simplest 
way to estimate Qm is to measure 
the skewness S3 of the convergence 
(Bernardea et al. 1997; Jain &; Sel- 
jak 1997; Schneider et al. 1998). We 



found that S3 was sensitive to the 
tails of the probability distribution 
function (pdf) and therefore required 
large sample sizes. Instead, we fit the 
pdf to an Edgeworth expansion and 
estimated S3 as a parameter. The re- 
sult is shown in figure M and shows 
that cosmological models with dif- 
ferent values of Qrn can be distin- 
guished at a high level of significance 
(the error bars are l-cr). The dashed 
curves show the predictions of per- 
turbation theory. 

Finally, figure ^ shows the con- 
vergence field centered on typical 
rich clusters in the Einstein-de Sitter 
model, rCDM, and the open CDM 
model. The fields are 10' arcminutes 
on a side and are chosen because cur- 
rent observational data is already ca- 
pable of detecting signal in the outer 
parts of clusters. Comparison of the 
EdS and open models shows that the 
convergence fields are morphologi- 
cally different: clusters in the open 
model appear more compact, regular 
and isolated. The clusters in the EdS 
model are not fully relaxed and ap- 
pear linked to the surrounding large- 
scale structure through filaments or 
more irregular structures. These dif- 
ferences can be quantified by using 
topological measures and by measur- 
ing moments such as the quadrupole 
or higher moments. In figure H we 
simply show the mean profile of the 
convergence around clusters, which 
also shows differences between the 
two models. The source galaxies have 
been taken to be at z = 1 and the 10 
richest clusters in a field 3 degrees on 
a side have been used to obtain the 
mean profiles. 
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Figure 2. The reconstructed convergence field from noisy ellipticity data on 
a single 3° field. The upper left panel shows the field without noise or sparse 
sampling. The upper right panel shows the reconstructed field using 2 x 10^ 
galaxies per square degree at a mean z = 1. The lower left panel used 1/Ath 
as many galaxies. The lower right panel shows a map of the "curl" field gen- 
erated from the same ellipticity data as in the lower left panel. Fluctuations 
in this field are solely due to noise, therefore comparison of the two panels 
shows the significance of the reconstructed features in the lower left panel. 



3 Conclusion 

We have shown results on estimat- 
ing the power spectrum and VLm from 
simulated, noisy weak lensing data. 
With several fields a degree on a side, 
the dark matter power spectrum on 
scales of 1 — 10 /i~^Mpc can be 
probed and Q^rn estimated to within 
about 0.1-0.2. We have also shown 
preliminary results on the morpho- 
logical differences between clusters in 



different cosmologies. Further work 
is needed to quantify the differences 
and explore the effects of noise and 
of varying the redshift distribution of 
source galaxies. 
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Figure 3. Convergence power spectrum estimated from simulated noisy data. 
The shear field on a single field 3° on a side is sampled by randomly dis- 
tributed galaxies with intrinsic ellipticites assigned as in figure ^. The power 
spectrum of the reconstructed n from the ellipticity data is shown with error 
bars obtained from 10 independent realizations of the noisy data. 



thank Anthony Banday for his pa- 
tient support in the writing of this 
contribution. 



References 

Bartelmann, M., Huss, H., Colberg, 
J.M., Jenkins, A. & Pearce, F.R. 
1998, A&A, 330, 1 

Bernardeau, P., van Waerbeke, L., &; 
MeUier, Y. 1997, A&A, 322, 1 

Blandford, R. D., Saust, A. B., 
Brainerd, T. G., &; Villumsen, J. V. 
1991, MNRAS, 251, 600 

Couchman, H.M.P., 
Barber. A. J. &: T homas. P.A. 1998, 
astro-p h/9810063| 

Gunn, J. E. 1967, ApJ, 147, 61 

Jain, B., & Seljak, U., 1997, ApJ, 
484, 560 Jain. B.. Seliak. U.. 



White, S. 1999, astro-ph/9901191 



Kaiser, N. 1992, ApJ, 388, 272 
Kaiser, N. 1998, ApJ, 498, 26 
Miralda-Escude, J. 1991, ApJ, 380, 1 



Premadi, P., Martel, H., Matzner, 
R., 1998, ApJ, 493, 10 

Schneider, P., Ehlers, J., & Falco, 
E.E., 1992, Gravitational Lensing 
(Springer Verlag, Berlin) 

Stebbins, A. 1996, |astro-ph/9609149| 

van Waerbeke, L., Bernardeau, P., 



MeUier, Y., 1998, ^stro-ph/9807007 

Wambsganss, J., Gen, R., Ostriker 
J. P. 1998, ApJ, 494, 29 



6 



Jain et al. 




-100 



10 

0[arcmin] 



Figure 4. The density parameter flrn estimated from the pdf of simulated 
noisy data. The skewness parameter S3 is estimated by minimizing the 
with respect to an Edgeworth expansion of the pdf. The four curves with 
decreasing peak heights are for the open, cosmological constant, Einstein- de 
Sitter models, all with V = 0.2 CD M power spectra, and an Einstein-de Sitter 
model with V = 0.5 CDM power spectrum. 



Numerical simulations of weak lensing 



7 





Figure 5. Clusters in open and Einstein-de Sitter cosmologies. The conver- 
gence in fields 10' on a side centered on a rich cluster is shown for the EdS 
model in the upper panels and for the open model in the lower panels. The 
values of the convergence range from over 10% in the center to below 1% in 
the outermost regions. Figure gives the mean profiles of the convergence. 
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Figure 6. The profile of the convergence in clusters in open and EdS cosmolo- 
gies. The dashed curve shows the average convergence profile measured from 
10 clusters in the EdS model while the solid curve shows the corresponding 
profile for the open model. 



